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Expansion by regions: 

revealing potential and Glauber regions automatically 
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Abstract When performing asymptotic expansions using 
the strategy of expansion by regions, it is a non-trivial task 
to find the relevant regions. The recently published Mathe- 
matica code asy . m automates this task, but it has not been 
able to detect potential regions in threshold expansions or 
Glauber regions. In this work we present an algorithm and 
its implementation in the update asy2 . m which also reveals 
potential and Glauber regions automatically. 



1 Introduction 

If a given Feynman integral depends on kinematic invariants 
and masses which essentially differ in scale, a natural idea 
is to expand it in ratios of small and large parameters. As a 
result, the integral is written as a series of simpler quanti- 
ties than the original integral itself and it can be substituted 
by a sufficiently large number of terms of such an expan- 
sion. For limits typical of Euclidean space (for example, the 
off-shell large-momentum limit or the large-mass limit), one 
can write down the corresponding asymptotic expansion in 
terms of a sum over certain subgraphs of a given graph [ 1- 
6]. This prescription of expansion by subgraphs has been 
mathematically proven (see [4] and Appendix B.2 of [6]). 
Moreover, there is an automated tool [7, 8] where such an 
expansion by subgraphs is implemented. 

For limits typical of Minkowski space (i.e. which cannot 
be formulated in Euclidean space) the universal strategy of 
expansion by regions [9-11, 6] is available. It consists of the 
following prescriptions: 

- Divide the space of the loop momenta into various re- 
gions and, in every region, expand the integrand in a 
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Taylor series with respect to the parameters that are con- 
sidered small there. 

- Integrate the integrand, expanded in the appropriate way 
in every region, over the whole integration domain of the 
loop momenta. 

- Set to zero any scaleless integral. 

As shown in [11], this prescription can also be applied to 
parametric representations of Feynman integrals, i.e. alpha 
parameters (or generalized Feynman parameters) integrated 
from to °°, eventually restricted by a delta function. Then 
the regions are specified by scaling relations between the 
parameters. 

There is no mathematical proof that this prescription is 
correct in all situations. But also no examples are known 
where a proper application of the expansion by regions leads 
to wrong results. An indirect proof exists for limits typical 
of Euclidean space because here the strategy of regions is 
equivalent to the mathematically proven expansion by sub- 
graphs. A systematic study of the expansion by regions was 
presented recently in [12]. There it was shown explicitly 
and illustrated using various one-loop examples that one 
can start from a decomposition of a given integral into non- 
intersecting domains and arrive at an expansion by regions 
in the above sense. This requires certain conditions on the 
choice and completeness of the considered regions which 
were derived in [12]. As pointed out there, the appearance 
of additional overlap contributions can be avoided by ade- 
quate choices of the regions and regularization parameters . 

While these findings provide some hints on the proper 
choice of the regions, it remains a non-trivial task to actu- 
ally reveal the typical regions for a given limit. Usually, one 
starts from considering one-loop examples, checks the re- 
sults against known analytical results, then proceeds in two 
loops etc. One can also use the second version [13] of the 
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code FIESTA [14] to obtain numerically several first terms 
of a given asymptotic expansion. 

Recently an algorithm for an automatic search of regions 
was suggested and implemented on a computer as the open 
source Mathematicacode asy .m [15]. The algorithm uses 
a geometric approach based on finding the convex hull of a 
set of points determined from a parametric representation of 
the Feynman integral. In this way all possible sets of scalings 
for the (Feynman) parameters are found which lead to non- 
vanishing (because non-scaleless) integrals. These regions 
may then be used to expand the parametric integral, or they 
can be translated into regions for expanding the integral in 
loop -momentum space. 

This code works successfully for a large variety of lim- 
its, at least in cases where the function J 7 in the correspond- 
ing parametric representation which depends on kinematic 
invariants and masses is only composed of terms with the 
same sign. Moreover, it was shown in [15] that in this case 
there are no regions except for the ones produced by the 
code (all others result in scaleless integrals). In particular, 
the code asy . m works for Sudakov-type limits which are 
typical of Minkowski space. As it was pointed out in [15] the 
code does not reveal potential regions in threshold expan- 
sions, and similarly it fails to detect the so-called Glauber 
regions. 

Our goal in the present paper is to provide an algo- 
rithm and the corresponding update asy2.m of the code 
asy . m which automatically identifies all regions relevant for 
a given integral, including potential and Glauber regions. We 
start in Section 2 by introducing the parametric representa- 
tions of loop integrals which we use later and by explaining 
how the contribution of a given region is obtained in the lan- 
guage of such parametric integrals. Then we elaborate our 
algorithm for revealing potential and Glauber regions and 
explain how the code asy2 .m is applied in such cases. This 
is done in Section 3 for an example with a potential region, 
and Section 4 deals with the more complicated problem of 
revealing Glauber regions. In both Sections 3 and 4, we first 
formulate simple changes of variables and decompositions 
of a given Feynman integral, using instructive one-loop ex- 
amples, which lead to integrals where asy2 . m is able to de- 
tect the relevant regions and print their list in terms of the 
scalings of the parameters. Then, for both cases, we explain 
how to use the new features of asy2 . m to perform these al- 
gorithmic steps automatically. 

In the case of Section 4 with Glauber regions, the struc- 
ture of the regions differs depending on whether the expan- 
sion is performed in loop-momentum space or at the level of 
the parametric integral. We show in Section 5 how to disen- 
tangle and match the various regions arising in this problem 
by using generic propagator powers, and how asy2 . m can 
be employed to automate such an analysis. 



A summary of the new features and the syntax of 
asy2 . m (together with a download link) is provided in Sec- 
tion 6. In Section 7 we conclude by discussing the mathe- 
matical problem of proving the expansion by regions for a 
simple example which is not related to Feynman integrals, 
but where asy2 .m works successfully. 

2 Expansion by regions in parametric representations 

We are dealing with dimensionally regularized Feynman in- 
tegrals 



F(q u ...,q„;a u ...,a N ;d) = 



1 



n^r, 

,=1 111=1^1 



, (1) 



where h is the number of loops, the indices a/ are general 
powers of the propagators, the dimension is d = 4 — 2e and 
the denominators E[ are given by 

I, I, 
Et = £ A) 3 k f -kj + ^B)- k, +D, + iO , (2) 

i.e. they are quadratic or linear functions of the external mo- 
menta qj and the loop momenta kj with the usual infinites- 
imal imaginary part +/0. Monomials in the numerator are 
taken into account as denominators raised to negative pow- 
ers. 

The alpha representation of (1) takes the form 

e -i1t{a+hd /2)/2 
F(q u ...,q n ;a u ...,a N ;d) = (i7t d/2 ) h N 

x jf°- • • J^f\ (da, a';i- l )u- d ' 2 e- iT l u , (3) 

where a = a\. The functions U and T depend polynomi- 
ally on the alpha parameters 0£/. Furthermore, U and J 7 are 
homogeneous functions of the alpha parameters with the ho- 
mogeneity degrees h and h + 1 , respectively. The function T 
is linear in the kinematic invariants and/or squared masses 
which we denote by S(, i = 1,2, . . ., while the function U is 
independent of the s;. 

If (1) is an integral with standard propagators 1 / (pf — 
mj + z'O) associated with the lines of a graph, then the func- 
tions U and T are called Symanzik polynomials and are 
given by the well-known formulae in terms of trees and 2- 
trees. For a general Feynman integral of the form ( 1) one 
can obtain these functions using the simple public code 
UF.m 1 [16] which is also part of the codes asy.m and 
asy2 .m. 



'The function UF[] from UF.m [16] is called with three arguments: 
The list of loop momenta, the list of denominators of the propagators 
and a list of replacement rules for all kinematic invariants. The output 
is a list with the following entries: the function U, the function T and 
the number of loops. In order to obtain U and J 7 with the correct sign, 
denominators have to be specified with the opposite sign as in (2), i.e. 
corresponding to a negative imaginary part — ;0: UF[{kl,k2, . . .}, 
{-El , -E2 ,...}, {replacement rales}] . 
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If some of the indices a\ are negative integers, i.e. they 
correspond to numerators instead of denominators of the in- 
tegral (1), the alpha representation (3) is to be understood in 
the limit where these indices tend to their negative integer 
values. Effectively, the integration over the corresponding 
parameters a/ is replaced by differentiating with respect to 
these parameters and setting them to zero. 

Closely related to (3) is the (generalized) Feynman para- 
metric representation 



F(qi,...,q n ;ai,...,a N ;d) = (m d > 2 ) 



_^ d /2, h e- i7ta r(a-hd/2) 

nf=ir(a ; ) 



dxi--- / dx N S\ Yxi-1 )I(xi,...,xn;si,S2,...) 
Jo Jo \f^ v I 



(4) 



where v in the argument of the delta function is an arbitrary 
non-empty subset of {1, . . . ,N}, 



I(xi, 



,x;v;si,S2, 



)=IP 



(5) 



and the functions U and J 7 are the same as those in (3) with 
the parameters 0£/ replaced by x/. It is well known that the 
formula (4) holds for any choice of the subset v in the ar- 
gument of the delta function. 2 This feature is related to the 
above-mentioned homogeneity properties of the functions U 
and J 7 . 3 If one chooses v = { 1 , . . . ,N}, the standard Feyn- 
man parametrization is recovered. 

Let us suppose that we have to study the asymptotic be- 
haviour in a one-scale limit, i.e. every mass and kinematic 
invariant has a certain scaling s, — > s' t = p K 'Sj, i = 1 , 2, . . ., ex- 
pressed in powers of the small parameter of the problem, p . 
The strategy of expansion by regions formulated in terms of 
parametric integrals (3) or (4) [11, 6] states that the asymp- 
totic expansion in such a limit is given by a sum over re- 
gions which are specified by the scalings of the parameters 
(Xi or xi expressed in powers r\ of the expansion parame- 
ter p . Each region r is labelled by the list r = {r\ , . . . , r^} 
of its scaling powers. The contribution of the region r is ob- 
tained by scaling the masses and kinematic invariants ac- 
cording to the given limit as specified above, by substitut- 
ing a/ -> a\ = p n ai , l = l...,N,in the integrand of (3) or 
x/ — > x'i = p n xi in the integrand of (4) and by expanding the 
integrand in powers of p . Here the product of the differen- 
tials da/ or dx{ provides another factor p^ in to the power 
counting. 



See e.g. the discussion in Section 3.4 of [17]. 
3 See also (36), (37) and Footnote 5 (p. 12) for a general proof. 



Explicitly, the contribution of the region r is given by the 
prefactorin (4) times p^ in times the integral 



dxi- - dx N 8 TV/-! )l(x' u ...,x' N ;s' u s' 2 ,...) (6) 

J \ftv J 



with the integrand expanded in powers of p. This expan- 
sion also involves the argument of the delta function in (6), 
such that, under the expansion, certain parameters drop out 
of the argument of the delta function and are integrated from 
to °°. For this reason the upper integration limit of all Feyn- 
man parameters should be kept at infinity and not switched 
to 1 even if, before the expansion, their integration is re- 
stricted by the delta function. One may avoid expanding the 
delta function by choosing the original subset v in (4) suffi- 
ciently small. 

Let us write down the leading-order (LO) contribution 
of a given region in a more explicit way. For the two basic 
functions in (5) we have 

''max 

U(x' 1) ...,x' N )= ]P p J Uj{x u ...,x N ) , 

J H min 
/max 

J r {x\,...,x' N ;s\,s' 2 ,...)= p j J r j(xi,...,x N ;s 1 ,s 2 ,...), 

j fm'm 

(7) 

where the arguments of the polynomials U and T on the 
left-hand side indicate that they are expressed in terms of 
the scaled parameters xj and s\, while their expansion coef- 
ficients on the right-hand side are expressed in terms of x; 
and Sj. According to the prescription formulated above, the 
LO contribution of the region r is represented as 



P 1/ n«l +'<mm («- (A+l) J ) +fmm(h 4 ~a) 

x{ . ndl2)h e-^r{a-hd,2) 



W =l r(a,) Jo 



1=1 v J 




In principle, the argument of the delta function in (8) con- 
tains a sum over only those scaled parameters x 1 , = p r 'x\ with 
the minimal scaling power r\ = r m ; n = minjri , . . . , rjv}. But 
by rescaling x/ — > p _rmin x/V/, the delta function is trans- 
formed into its standard form without powers of p, while 
the rest of the integral remains invariant due to the homo- 
geneity of the polynomials U Umin and Ff min ■ Finally, as for 
the original Feynman parametric representation (4), one can 
choose again an arbitrary non-empty subset Vo for the sum 
in the argument of the delta function in (8). 

The list of scalings r = {n , . . . , r^} of a region is only 
determined up to adding the same arbitrary real number c to 
each entry, because the corresponding contribution stays the 
same under r/ — s- r\ +cV/. In particular, the LO behaviour 
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presented in (8) is independent of such a shift c, because 
Mmin -> "min + he and / min -> / m i n + (/i + 1 )c, due to the ho- 
mogeneity properties of U and T. 

The terms of the expansion come from various regions 
and can be ordered according to accompanying powers of p. 
After keeping some first terms of the expansion one can set 
p = 1 and write down the given Feynman integral as these 
selected first terms plus a remainder which vanishes suffi- 
ciently fast in the given limit. 

It turns out that only a limited number of regions con- 
tribute to the expansion because for the rest of the regions 
one obtains integrals without scale which are set to zero. It 
is the subject of this paper and the task of the code asy . m 
and its updated version asy2 . m to find all relevant regions 
for a given integral. 

3 Revealing potential contributions 

Let us consider the one-loop propagator diagram with two 
massive lines in the threshold limit, i.e. when y = m 2 — 
q 2 /4 — > with q being the external momentum: 

F(q 2 ,m 2 ) = 



d d k 



(k 2 -m 2 ) {(k-q) 2 



(9) 



where the usual +/0 is implied in all the propagators. Within 
the strategy of expansion by regions, the hard and the poten- 
tial regions give contributions to the expansion [9, 6]. The 
previous version of the code asy .m [15] reported only about 
the hard region. The reason for this can be seen in the corre- 
sponding parametric representation, 

F(q 2 ,y) = in d / 2 r(e) 



(x\ + x 2 ) 8 (x\ +X2 — 1) dxidjQ 



i-(x l -x 2 ) 2 +y(x 1 +x 2 ) 2 -i0 



(10) 



where the parameters x, are integrated from to °° (restricted 
by the delta function). As it was pointed out in [15], it is 
the region where x\ w x 2 (more precisely x\ — x 2 ^y 1 ^ 2 ) 
which causes problems. In other words, the polynomial in 
the square brackets in (10) (considered at positive q 2 and 
y) has terms of different sign, such that cancellations occur 
because of the presence of the negative term — q 2 x\X2/2. 

To reveal the missing potential contribution, let us per- 
form a simple trick. We decompose the integration domain 
into two subdomains, x\ < X2 and X2 <X\. The two resulting 
integrals are equal to each other, but such an equality will 
not generally take place for any integral. In the first domain 
we turn to new variables by x\ = x J l /2, x 2 = x' 9 + x\/2, re- 
move the primes at x; and obtain the integral (again from 
to °° with the usual restrictions via the delta function) 

r(e) ff (xi +x 2 ) 2e ~ 2 8{xi + x 2 - 1) dxidx 2 



lit 



1/2 ■ 



4 A 2 



+y(x l +x 2 ) 2 -i0 



(11) 



The goal of this trick was to make the line x\ = X2 (in the old 
variables) the border of an integration domain which turned 
out to be (in the new variables) X2 = 0. Now we can ran 
the code asy2 . m. Since this is a parametrical integral rather 
than a Feynman integral we use the newly introduced com- 
mand 4 WilsonExpand [] for integrals where all parameters 
are integrated from to °°: 

WilsonExpand [q~2/4*x2~2 + y*(xl + x2)~2, 
xl + x2, {xl, x2>, {q -> 1, y -> x} , 
Delta -> True] 

The first two arguments of WilsonExpand [] are the poly- 
nomials J 7 and U, respectively, as defined in Section 2. They 
can easily be determined from the square brackets in the 
denominator of the parametric integral (11) and from the 
round brackets in the numerator. The third argument is the 
list of integration parameters, and the fourth argument speci- 
fies the scaling of the kinematic quantities with respect to the 
small parameter which is labelled by the global symbol x. 
Here by y -> x we tell the code that y is the small expan- 
sion parameter, and by q -> 1 we specify that the momen- 
tum q scales as y° = 1. The option Delta -> True tells 
WilsonExpand [] that, under the integration, the sum over 
an arbitrary non-empty subset of the integration parameters 
is restricted to 1 by a delta function. 

Note that WilsonExpand [] can only take into account 
such a delta function if the specific choice of the sum over 
parameters in the argument of the delta function is irrelevant. 
This is the case for the generalized Feynman parametric in- 
tegral (4) introduced in Section 2. The integrals (10) and (1 1) 
are special cases of (4) such that we could e.g. replace 
8(x\ +X2 — 1 ) by 8(x\ — 1 ) without changing the value of the 
integrals. If, however, a specific form of the delta function is 
assumed, e.g. by replacing x\ +X2 —> 1 under the integral, 
then the option Delta of WilsonExpand [] does not apply 
(see Section 6 for details). Alternatively, WilsonExpand [] 
can be used without the option Delta after eliminating one 
of the integrations, e.g. via 8(x\ — 1). 

The output of the above-mentioned call of 
WilsonExpand [] is a list of regions specified by the 
scaling of the parameters x, in powers of the small 
parameter y: 

{{0, 0}, {0, 1/2}} 

This specification of a region corresponds to the list of scal- 
ings {rj , . . . , r^} introduced in Section 2. The first entry in 
the output, {0,0}, refers to the hard region with the scaling 
xi ~ y°, X2 ~ y° which has already been found by asy .m. 
But now also the region {0, 1/2} is found with the scaling 
xi ~ y°, X2 ~ y x l 2 which provides the potential contribution. 
The contribution of the hard region starts with the order 



y 



Every term of the expansion can be evaluated in terms of 



4 The name of the command refers to its application to parametric inte- 
grals contributing to Wilson loops. 
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gamma functions for general e. According to the prescrip- 
tions for writing down the contribution of a region formu- 
lated in Section 2, the contribution of the k-th order expan- 
sion of (1 1) in the potential region reads 



lit 



d/2 



2k\ 



ax\ 



d \ k x* e - z 8(xi-l) 



Only the leading order (k = 0), 



lit 



ill! 



dx 2 



t4+y) 



(12) 



(13) 



yields a non-vanishing contribution which is evaluated in 
terms of gamma functions at general e. Taking into account 
that we have two identical integrals after our decomposition, 
we arrive at the following result for the potential contribu- 
tion which is of order 

y l/Z-e. 



in d ' 2 r{t 



(14) 



in agreement with [6]. 

In fact, such a trick of making manifest squares of some 
linear combination of the integration parameters was already 
used in the code FIESTA [14] in order to evaluate numer- 
ically Feynman integrals at a threshold. Using the imple- 
mentation of this procedure in FIESTA it turned out to be 
possible to automate the above trick for a general Feynman 
integral. In the present version asy2.m, the user may call 
the command AlphaRepExpand [] with the additional op- 
tion PreResolve enabled which automatically looks for the 
change of variables described above: 

AlphaRepExpand [{k} , 

{k~2 - m~2, (k - q)~2 - m~2}, 
-Cq-2 -> qq, m~2 -> qq/4 + y>, 
-Cqq -> 1, y -> x}, 
PreResolve -> True] 

As in the previous version of asy.m, the arguments of 
AlphaRepExpand [] are the list of loop momenta, the list 
of denominators of the loop integral, a list of replacements 
for the kinematic quantities, and the list of scalings with re- 
spect to the small parameter x. The output is a list of entries 
which indicate for each region the changes of variables, the 
Jacobian of the integral transformation, and the scalings of 
the new variables: 

«{x[l] -> y[l] + y[2]/2, x [2] -> y[2]/2}, 

2, -CO, -1/2}}, 
«x[l] -> y[l] + y[2]/2, x[2] -> y[2]/2}, 

2, {0, 0}}, 
{{x[l] -> y[l]/2, x[2] -> y[l]/2 + y[2]}, 

2, {0, 0}}, 
«x[l] -> y[l]/2, x[2] -> y[l]/2 + y[2]}, 

2, -CO, 1/2}}} 



The original Feynman parameters are labelled x [1] , x [2] , 

. . ., while the new parameters are labelled y [1] , y [2] , 

The scaling relations which determine the regions are spec- 
ified for the new parameters. 

As explained in Section 2, a region remains invariant 
if all of its scalings are shifted by the same amount, i.e. if 
all parameters are rescaled simultaneously. So the scaling 
{0,-1/2} shown above is equivalent to {1/2,0}, i.e. the 
first new parameter is suppressed by y 1 / 2 with respect to the 
second one. For both parts of the integral decomposition the 
hard and potential regions are found. 

The preresolution algorithm implemented in asy2.m 
(and switched on with the option PreResolve) tries to elim- 
inate factorized combinations of terms in the function T 
which potentially cancel each other, like {x\ ~ x 2 ) 2 in the 
example above. It checks all pairs of variables (say, x and 
y) which are part of monomials with opposite sign. For all 
those pairs the code tries to build a linear combination z of 
x and y such that in the variables x and z or y and z this 
monomial disappears. The code checks whether in the new 
variables the number of monomials with opposite sign de- 
creases. For all such pairs the code recursively repeats the 
initial procedure in the new variables. As a result it creates a 
tree of possible bisections and corresponding replacements 
of variables. A leaf of this tree is a set of sectors and func- 
tions such that one cannot decrease the number of monomi- 
als with opposite sign any longer. Ideally it means that all 
monomials now have the same sign. The code analyzes all 
leafs and chooses one of those with the minimal number of 
opposite-sign monomials (or the minimal number of sectors 
if the numbers of monomials with opposite sign coincide). 
After finishing with the preresolution, the code performs the 
replacements and looks for regions in all those sectors, using 
the algorithm of the original code asy . m described in [15]. 

Note that the algorithm can only find the necessary vari- 
able transformations if it is able to determine the relative 
signs of all terms in the polynomial T. As the signs of sym- 
bols are unknown to asy2 . m, the substitution rules of the 
third and fourth arguments of AlphaRepExpand [] must re- 
place all kinematic quantities by numbers (integers or frac- 
tions of integers) or powers of the small parameter x. 

We have checked that the updated version asy2 . m works 
in various examples of the threshold expansion (considered 
in [9, 6]): a triangle, a box, the two-loop propagator dia- 
gram (with the masses m,m,m,m,Q), a two-loop vertex di- 
agram. Because of the decomposition of a given integration 
domain into subdomains, the number of resulting integrals 
for various regions increases a little bit. For example, the 
(hard-hard) region for the two-loop propagator diagram is 
described by six integrals, the (potential-ultrasoft) region is 
also described by six integrals, etc. However, the (potential- 
hard) region is described by four integrals with some regions 
(with scalings composed of powers 1, 1/2 and 0), and four 
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Pi +P2 



9i +92 




Fig. 1 One-loop five-point integral exhibiting a Glauber contribution 

more integrals with a set of regions of a different type (com- 
posed of 1 and 0). 

Let us finally mention that the preresolution algorithm 
also works for threshold expansions with unequal masses. 
Returning to the one-loop example (9), but now with differ- 
ent masses m\ and mi in the propagators and the expansion 
parameter defined as y = (mi + m2) 2 /4 — q 2 /4- — s- 0, we call: 

AlphaRepExpand [{k} , 

{k~2 - ml~2, (k - q)~2 - m2~2}, 
{q~2 -> (ml + m2)~2 - 4*y>, 
-Cy -> x, ml -> 2, m2 -> 3/2}, 
PreResolve -> True] 

We have set the values of the masses to different ratio- 
nal numbers, mi 2 and mi —> 3/2 (the actual values are 
irrelevant), thus permitting the preresolution algorithm to 
know about their sign without assuming any equality or 
other relation between them. The polynomial T in the para- 
metric representation of the Feynman integral reads T = 
(mixi — mixi) 2 + Ayxixi, and cancellations occur where 
mixi = mixi- These cancellations are automatically made 
explicit by adequate changes of variables. The output 

{{{x[l] -> y[l] + 3*y[2]/7, 

x[2] -> 4*y[2]/7}, 7/4, {0, -1/2}}, 

{{x[l] -> y[l] + 3*y[2]/7, 

x[2] -> 4*y[2]/7}, 7/4, {0, 0}}, 

{{x[l] -> 3*y[l]/7, 

x[2] -> 4*y[l]/7 + y[2]}, 7/3, {0, 0}}, 

{{x[l] -> 3*y[l]/7, 

x[2] -> 4*y[l]/7 + y[2]}, 7/3, {0, 1/2}}} 

shows that asy2 . m detects both regions. The variable trans- 
formations are always normalized such that sums of the pa- 
rameters (here xi +X2) remain invariant. 

So, we now have a manifestly Lorentz-invariant treat- 
ment of threshold expansion and a code that automatically 
provides the set of relevant regions. 

4 Revealing Glauber contributions 

Let us consider the one-loop five-point integral in Fig. 1, 
where two initial-state partons both perform a collinear split- 
ting into two partons each with momenta pi,p2 and qi,q2, 



respectively. While two partons, one of each pair, collide 
with a large centre-of-mass energy Q = y/ (p2 + qi) 2 , the 
two remaining partons exchange a particle with the small 
mass m. We will use the simplified kinematics p\ = p2= p 
and q\= q2 = q with p 2 = q 2 = and (p + q) 2 = 2p ■ q = Q 2 
in the limit m 2 /Q 2 0: 



F(Q\m 2 ) = 



A d k 



(k 2 — m 2 )(k 2 — 2p ■ k)(k 2 + 2p ■ k) 
1 



(k 2 ~2q-k)(k 2 + 2q-k) 



(15) 



Before we search for regions using asy . m, we notice 
that this five-point integral is similar to the Sudakov form 
factor example treated in Section 6 of [12]. From the view- 
point of the convergence of the expansions, the second and 
third propagators of ( 15) are equivalent, and so are the fourth 
and fifth propagators. Effectively, the five-point integral has 
only three different types of propagators which are equiv- 
alent to the ones of the three-point integral in [12]. So the 
integral (15) can be expanded in loop-momentum space 
employing the regions and convergence domains known 
from [12] (and using generic propagator powers as analytic 
regulators where necessary): 

- a hard region where k ~ Q, 

- a 1 -collinear region where k 2 ~ p-k^ m 2 and q-k~ Q 2 , 

- a 2-collinear region where k 2 ~ q-k^m 2 and p-k~ Q 2 , 

- a Glauber region where p ■ k ~ q ■ k ~ m 2 , and the com- 
ponents of k perpendicular to the plane spanned by p,q 
scale as k± ~ m. 

The collinear-plane region mentioned in [12] yields only 
scaleless contributions. But, in contrast to the three-point in- 
tegral, the five-point integral has a non-vanishing Glauber 
contribution. The Glauber region even provides the lead- 
ing contribution scaling as (m 2 )~ 2 ~ £ , whereas the collinear 
contributions start with (m 2 ) _1_£ and the hard contribution 
starts with (f7z 2 ) . 

The five-point integral (15) can be represented in terms 
of an integral over Feynman parameters, 



F(Q 2 ,m 2 



-;V //2 r(3 + e) /••• /djci-"djC5 
8(Z, Xi -l) (x 1 + ...+x5) 1+2£ 



[xi(xi + ... +X5)m 2 + (x2 — xt,)(xa, — xs)Q 2 — ;0 



3+£ 



(16) 



where one can choose the sum in the argument of the delta 
function in an appropriate way, i.e. restrict only the sum over 
a subset of the parameters to 1 and extend the integration 
over the rest of the parameters to the whole domain [0,°°). 

Applying the strategy of expansion by regions in 
Feynman-parameter space and trying to reveal regions rele- 
vant to the given limit m 2 /Q 2 — s- with the help of the code 
asy .m [15], we call: 
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AlphaRepExpand [{k} , 

{k~2 - m~2, k~2 - 2*p*k, k~2 + 2*p*k, 

k~2 - 2*q*k, k~2 + 2*q*k}, 
-Cp-2 -> 0, q~2 -> 0, p*q -> Q-2/2}, 
{Q -> 1, m~2 -> x}] 

The output states the following set of three regions: 

{{0, 0, 0, 0, 0}, {0, 0, 0, 1, 1}, 
{0, 1, 1, 0, 0» 

As before, the regions are specified by the scaling of the 
Feynman parameters in terms of powers of the small pa- 
rameter m 2 . For example, for the second region we have 
x\ ~ X2 ~ xt, ~ (in 2 ) , X4 ~ X5 ~ (in 2 ) 1 ■ The first region 
is hard; its contribution starts with (in 2 ) . The second and 
third regions start with order (m 2 )~ l ~ £ . They correspond 
to the two collinear regions stated for the momentum-space 
expansion above. But asy .m does not find anything corre- 
sponding to the Glauber region; in particular, none of the 
regions found by asy . m provides the leading (m 2 )~ 2 ~ £ con- 
tribution. 

We notice that, as in the previous section about po- 
tential contributions, the polynomial in the square brackets 
of (16) has terms of different sign. The missing Glauber 
contribution stems from the parameter region where either 
(x2 — X3) ~ (m 2 ) 1 or (x4 — X5) ~ (m 2 ) 1 . So let us decompose 
the parametric integral into four parts corresponding to the 
domains where the two factors (x2 — X3) and (X4—X5) are ei- 
ther positive or negative and then introduce new variables in 
such a way that this product takes the form ±x 2 x / 4 . For exam- 
ple, in the domain X2 < X3, X5 < X4 we change the variables 
by X2 = X3/2, X3 = x' 2 +x' 3 /2 and by X4 = x' 4 +x' 5 /2, x$ = 
X5/2, similarly to our example in the previous section. How- 
ever, in the threshold expansion the cancelling terms ap- 
peared in squared form such that a transformation between 
one pair of variables was sufficient. Here two separate fac- 
tors involve cancellations, which requires a twofold change 
of variables. 

Removing the primes from the variables xi, the paramet- 
ric integral reads 

F(Q 2 ,m 2 )=2(I + +I_) (17) 
with 

I ± = -i^ r -^J...j6x^ 

5 (jfl — 1 ) {x\ +X2+X3 +X4+X 5 ) 1+2E 

[x\ (x\ +X2 +X3 + X4 + xs)m 2 ± X2X4Q 2 — ;0] 3+£ 

(18) 



where we have chosen the argument of the delta function as 
x\ — 1, so that we may also write 

/±= _^/2£^±£) r... r<ix 2 -dx 5 

4 Jo Jo z 

(1 +x 2 +x 3 +x 4 +x 5 ) 1+2e 
x i - ^ 3— . (19) 

[(1 +X2 +X3 +X4 +x$)m 2 ±X2X4Q 2 — iQ]~ 

It is sufficient to consider the expansion of /+ and obtain 
a result for /_ by analytically continuing Q 2 — > — Q 2 — /0, 
taking into account that the dependence on Q 2 is power-like. 

Now we can apply asy2 . m to the integral I + using either 
the command 

WilsonExpand[xl*(xl + x2 + x3 + x4 + x5)*m~2 
+ x2*x4*Q~2, 
xl + x2 + x3 + x4 + x5, 
-Cxi, x2, x3, x4, x5>, 
{Q~2 -> 1, m~2 -> x}, Delta -> True] 

for integrals (18) restricted by a delta function (see Sec- 
tion 6), yielding the output: 

{{0, 0, 0, 0, 0}, {0, 1, 0, 0, 0}, 
{0, 0, 0, 1, 0» 

Or we use the command 

WilsonExpand[(l + x2 + x3 + x4 + x5)*m~2 
+ x2*x4*Q~2, 
1 + x2 + x3 + x4 + x5, 
{x2, x3, x4, x5} , 
-CQ-2 -> 1, m~2 -> x}] 

for integrals (19) over variables from to °° without any 
restriction, and obtain the output: 

{{0, 0, 0, 0}, {1, 0, 0, 0}, -CO, 0, 1, 0» 

The two results are equivalent, as x\ ~ (m 2 ) is implied in 
the second case. So we obtain again three regions. We will 
see in a moment that this list of regions is indeed correct and 
complete. 

But first, let us emphasize that the new preresolution al- 
gorithm in asy2 .m is capable of performing the transforma- 
tion of the integral from (16) to (18) automatically: 

AlphaRepExpand [{k} , 

{k~2 - m~2, k~2 - 2*p*k, k~2 + 2*p*k, 

k~2 - 2*q*k, k~2 + 2*q*k>, 
-Cp-2 -> 0, q~2 -> 0, p*q -> Q-2/2}, 
{Q -> 1, m~2 -> x}, 
PreResolve -> True] 

The output of this command lists four different variable 
transformations according to the twofold decomposition de- 
scribed above. For each of the integrals over new parame- 
ters, the regions {0, 0, 0, 0, 0}, {0, 1 , 0, 0, 0} and {0, 0, 0, 1 , 0} 
are found, up to permutations in the order of the parameters 
from different changes of variables. 



X 



The evaluation of the contributions to each region, as 
found by WilsonExpand [] or AlphaRepExpand [] (in- 
cluding the PreResolve option), is straightforward. The 
first region is the hard one. The contributions of the second 
and third regions are not individually regularized by dimen- 
sional regularization, as it often happens for Sudakov-type 
limits. We use an auxiliary analytic regularization by intro- 
ducing additional powers x^x^x^x^ 5 of the new variables 
into the integrand of (19), taking the limit 82, 83, 84, 85 — >• 
in the end. The leading-order (LO) contribution of the sec- 
ond and third regions to the integral F(Q 2 ,m 2 ) reads 



-IK 



d/2 _ 



inr{e) 

2Q 2 (m 2 ) 2 + £ 



(20) 



This agrees with the leading contribution of the Glauber re- 
gion in the momentum-space expansion. 

We have found the leading Glauber contribution of or- 
der (m 2 )~ 2 ~ £ . But we seem to have lost the two collinear re- 
gions with the scalings {0,0,0,1,1} and {0,1,1,0,0} found 
before the change of variables. In fact, we can evaluate the 
contributions from these two regions by expanding the inte- 
gral (19). The resulting integrals are scaleless and regular- 
ized by the parameters 83,85,50 they vanish, and asy2 . m is 
right in omitting these two regions. 

We are also able to solve the integral (19) including the 
auxiliary analytic regularization factor x^x^x^x^ in terms 
of a onefold Mellin-Barnes representation: 



-lit 



d/2 r(i + 8 3 )r(i + 8 5 ) 



-3-e- 



x—jdz(m'y(±Q z -M) 

x r(- z )r(-2 - e + 8 2 - z)r(-2 - e + 84 - z) 

r(l-82-8 3 -84-8 5 + z)r(3 + e + z) 
X T(-l-2e-z) 



(21) 



The asymptotic expansion of I± in the limit m 2 /Q 2 —> is 
obtained by taking the residues of the poles of the functions 
F(. . . — z). The poles of r(—z) correspond to the hard re- 
gion, while the poles of the two functions F (— 2 — e + 82 4 — 
z) provide the contributions of the second and third regions. 
So asy2 .m has found all contributing regions. 

In the Mellin-Barnes integral (21) we can safely take the 
limit 82,83,84,85 — > 0, add up I + and /_, and arrive at the 
Mellin-Barnes representation 

F(Q 2 ,m 2 ) = in*! 2 - — [ dz(m 2 Y(Q 2 )- 3 - £ -^ £ +^ 2 
2 2ni J 

,r(- z ,r(-2- e - z) r(^) r(>+;)r(^) 



r(-l-2e-z) 



(22) 

The LO contribution to F(Q 2 ,m 2 ) is obtained from the 
residue of the single pole at z = — 2 — e, in agreement 



with (20). The next-to-leading-order (NLO) contribution 
stems from the residue of the double pole at z = — 1 — £ and 
reads 

Q 2 



In- 



(23) 



This agrees with the NLO contributions of the second and 
third regions. 

At next-to-next-to-leading order (NNLO) there is a con- 
tribution from the residue of the single pole at z = — £ which 
reads 



-iTt 



1/2 



iitr(2 + e) 
4(Q 2 )\m 2 ) £ 



(24) 



and agrees with the NNLO contributions of the second and 
third regions. The second NNLO contribution comes from 
the residue of the single pole at z — 0. It is given by 



-lit 



d/2 



ie^/ 2 r(-2 - E )r r (ifs) 

2(g 2 )3+*r(-l-2e) 



(25) 



and agrees with the LO contribution of the hard region. 
So indeed all contributions to the five-point integral up to 
NNLO are correctly reproduced by the contributions of the 
three regions found by asy2 . m after the decomposition of 
the integral and the change of variables. 

When revealing Glauber regions for a general diagram, 
the preresolution algorithm of asy2 . m tries to eliminate 
monomials with opposite sign in the polynomial T by au- 
tomatically separating the integration into domains and per- 
forming changes of variables. If the option PreResolve is 
enabled for AlphaRepExpand [] , the code warns the user 
if the elimination of monomials with opposite sign has not 
been successful, such that possibly not all regions are re- 
vealed. This is the case if some monomials of opposite sign 
remain in the polynomial Rafter tries to eliminate them, or if 
symbols with unknown signs are present in the polynomial. 
We are therefore convinced that AlphaRepExpand [] , with 
the preresolution enabled, either reveals all relevant regions 
or issues a warning. 

As for the threshold expansion of the previous section, 
asy2 . m can also treat more complicated kinematical situ- 
ations, e.g. the five-point integral depicted in Fig. 1 with 
Pi 7^ P2 an d 1\ 7^ 12 (retaining p\ |j P2 and q\ \\ q2 such 
that pi ■ pj = q { qj = and 2p t ■ qj = P t Qj with P u Qj > 0): 

AlphaRepExpand [{k} , 

{k~2 - m~2, k~2 - 2*pl*k, k~2 + 2*p2*k, 

k~2 - 2*q2*k, k~2 + 2*ql*k}, 
{pl~2 -> 0, p2~2 -> 0, pl*p2 -> 0, 
ql~2 -> 0, q2~2 -> 0, ql*q2 -> 0, 
pl*ql -> Pl*Ql/2, pl*q2 -> Pl*Q2/2, 
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p2*ql -> P2*Ql/2, p2*q2 -> P2*Q2/2}, 

{m~2 -> x, PI -> 1, P2 -> 3, 

Ql -> 2, Q2 -> 3/2}, 
PreResolve -> True] 

All kinematic invariants are replaced by rational numbers 
in order to enable the preresolution algorithm to work. The 
code correctly decomposes the integral into four pieces and 
finds the three regions for each of them. 

Let us finally discuss the reason why, besides the hard 
region which is always present, the expansion in loop- 
momentum space requires two collinear regions and one 
Glauber region, whereas the expansion of the parametric in- 
tegrals (19) has two regions providing the leading Glauber 
contribution and no further collinear regions. In fact, the 
momentum-space expansion is also valid for loop inte- 
grals (15) where each propagator is raised to an arbitrary, 
even non-integer power. For the decomposition of the para- 
metric integral and the change of variables, however, we 
have assumed the specific form (15) with each propagator 
present exactly once. In this case, the loop integrand can be 
expanded into partial fractions as follows: 

1 

(k 2 - m 2 )(k 2 -2p-k)(k 2 + 2p ■ k) (k 2 -2q-k) (k 2 + 2q-k) 



1 



1 



1 



4(m 2 ) 2 \k 2 -m 2 k 2 (k 2 ) 2 , 
1 1 \ / 1 



k 2 -2p-k k 2 + 2p-kj \k 2 -2q-k k 2 + 2q-k 



(26) 



Expanding this product of terms, one obtains twelve three- 
point integrals, which are well known. Because they only 
depend on Q 2 = 2p ■ q, we recognize from the last two fac- 
tors the structure F(Q 2 ,m 2 ) = 2(1+ +I-), where 1+ and I 
are related by Q 2 — > — Q 2 as before. The three-point integrals 
with the massless propagators 1/k 2 or 1/ (k 2 ) 2 only have a 
hard region. The massive three-point integral with propaga- 
tor 1 /(k 2 — m 2 ) is known to possess a hard and two collinear 
regions. Its LO and NLO hard contributions are cancelled 
by the massless three-point integrals, such that the uncan- 
celled hard contributions start with (m 2 ) 2 / (m 2 ) 2 — (m 2 ) , 
as for the five-point integral. The LO collinear contributions 
of the three-point integrals, enhanced by the 1 /(m 2 ) 2 pre- 
factor, scale as (m 2 )~ 2 ~ £ . 

So, in the special case when all propagator powers are 
equal to 1, the five-point integral reduces to a linear combi- 
nation of three-point integrals revealing the same structure 
of regions as found from the expansion of the parametric in- 
tegral (19). This picture changes when generic propagator 
powers are introduced as analytic regulators, which is done 
in the next section. 



5 Disentangling regions via propagator powers 

In the previous section we have seen different patterns of 
regions arising when expanding either in loop-momentum 
space or in parametric space. The individual contributions 
can be disentangled more easily when the dependence on 
the propagator powers is retained. Instead of (15), let us con- 
sider the integral 



F(Q\m 2 ) 



d d k 



(k 2 -m 2 y+^(k 2 ~2 P -ky+ x i 



(k 2 +2p-k) (k 2 -2q-k) l +^ (k 2 + 2q-k) 



(27) 



where the analytic regularization parameters A, make the 
propagator powers different from the previous case. The 
asymptotic expansion in loop-momentum space yields con- 
tributions from the four regions listed in the beginning 
of Section 4. The LO hard contribution still scales as 
(m 2 ) , but now the LO 1 -collinear contribution scales 
as (ff^ 2 )~ I ~ £ ~^ 1 ~' 1 • 2 ~' , ' 3 , the LO 2-collinear contribution as 
(m 2 )~ l ~ e ~ Xl ~^~ Xs and the LO Glauber contribution as 
(m 2 )~ 2 ~ £ ~ Xl ~^ 2 ~^~ XA ~^ . So all regions are characterized 
by a distinct scaling. When performing the expansion in 
parametric space for generic A,-, we are able to disentangle 
the contributions from each region and match them with the 
regions obtained in loop-momentum space. 

For generic Ai,...,A5, the parametric integral corre- 
sponding to (16) reads 



F (Q 2 ,m 2 ) = - in d/ 2 e -i*{h+:.+h) 
r(3 + e + A 1 + ... + A 5 ) 



x \x\ (x\ + . . . +xs)m 



+ (x 2 -x 3 )(x 4 -x 5 )Q 2 -iO] 



2 ;nl -(3+e+^+...+%) 



(28) 



The change of variables performed in Section 4, e.g. x 2 = 
x' 2 +JC3/2, .£3 = JC3/2 forx2 > X3, is complicated by the pres- 
ence of the factors Xj 2 ^ 3 ^ 4 ^ 5 , where parts of the mono- 
mials will change into polynomials. If we want to keep the 
simple form F(Q 2 ,m 2 ) = 2(1+ +I-) from the previous sec- 
tion, then we have to require A3 = A2 and A5 = A4. Under 
this restriction of the parameters A,-, the parametric integrals 
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can be written as 



. dn e~ iK ^ +2h+2h) p(3 + e + Ai + 2X 2 + 2X4) 

±__? 4 l +^+hr(i + x l )r 2 {i + A 2 )r 2 (i +A 4 ) 



/ djci • • -dz 5 5 {xi - 1) x A Y l (2x 2 +x 3 ) X2 x\ 2 

(2X 4 +X 5 )*4j£ (JC1 + . . . +x 5 ) 1 + 2£+;i l+ 2;i 2+ 2;i 4 

3+£+A!+2A 2 +2A 4 



-xs)m 2 ±X2X4Q 2 — ;0] 



(29) 



In order to find the regions for the asymptotic expansion 
of (29), we have to provide the additional polynomial fac- 
tors (2x2 +^3) and (2x4 +X5) to WilsonExpand [] . We can 
do this by multiplying the new polynomials to the second 
argument of the command: 

WilsonExpand [xl*(xl + x2 + x3 + x4 + x5)*m~2 
+ x2*x4*Q~2, 
(xl + x2 + x3 + x4 + x5) * (2*x2 + x3) 

* (2*x4 + x5) , 
-Cxi, x2, x3, x4, x5>, 
{Q~2 -> 1, m~2 -> x}, Delta -> True] 

The output of this command reads: 

{{0, 0, 0, 1, 0>, {0, 0, 0, 1, 1}, 
{0, 1, 0, 0, 0}, {0, 1, 1, 0, 0}, 
{0, 0, 0, 0, 0» 

In addition to the regions present in the analysis of Sec- 
tion 4, we retrieve the two collinear regions with scalings 
{0,0,0,1,1} and {0,1,1,0,0}. The updated code asy2.m 
is capable of taking into account generic propagator pow- 
ers automatically when told through the additional option 
GenericPowers: 

AlphaRepExpand [{k} , 

{k~2 - m~2, k~2 - 2*p*k, k~2 + 2*p*k, 

k~2 - 2*q*k, k~2 + 2*q*k>, 
{p~2 -> 0, q~2 -> 0, p*q -> Q~2/2}, 
{Q -> 1, m~2 -> x}, 

PreResolve -> True, GenericPowers -> True] 

When the option GenericPowers is enabled, the polyno- 
mial U obtained from the loop integral is multiplied by 
the product of all Feynman parameters, x\X2 • • • , before the 
changes of variables are performed. Some of these addi- 
tional factors then turn into polynomials through replace- 
ments in the preresolution algorithm, while others remain 
monomials and are therefore irrelevant for the analysis of 
asy2 .m. The call of AlphaRepExpand [] stated above cor- 
rectly yields all five regions for each of the variable trans- 
formations. 

For the evaluation of the Glauber contributions with 
scalings {0,0,0,1,0} and {0,1,0,0,0}, an additional ana- 
lytic regularization is needed, and we choose to multiply the 



integrand of (29) by x^x^ 4 . (The parameters £3 , 5s from Sec- 
tion 4 are not needed here due to the presence of A2,A4.) 
The two Glauber contributions are individually singular in 
the limit 82,84 — > 0, but this singularity cancels in the sum 
of the two contributions. 

The leading contribution to the integral F(Q 2 ,m 2 ) origi- 
nates from the sum of the LO Glauber contributions. It reads 



lit 



ie -i7t(Xi+2h.+2Xn) r (1 



-A 2 )r(i + A 4 ) 



2Q2 (j n 2y+£+X l +21 2 +2X4 

r(2 + £ + Ai+2A 2 + 2A 4 )r(-l -£-2A 2 -2A 4 ) 

r(i +A0r(i + A 2 )r(i + A 4 )r(i - E ) 



(30) 



This agrees with the LO contribution from the one Glauber 
region in the momentum-space expansion of (27) for A3 = 
A2 and A5 = A 4 . Here, in the expansion of the parametric in- 
tegrals, we have two regions producing the Glauber contri- 
bution. This is possible because the contributions from both 

regions have the same scaling (for 82 = 84 = 0), starting 
with ( m 2 )- 2 -^-h-2X 2 -2l^ For A . = Qi the resuh ( 3 ) re _ 

produces (20). 

Among the NLO contributions to F(Q 2 7 m 2 ), we expect 
NLO Glauber contributions and LO collinear contributions. 
However, the NLO Glauber contributions vanish exactly for 
general X\, X2 = A3 and A4 = A5, due to non-trivial can- 
cellations between the pieces which contribute to the NLO 
expansion of (29) for either of the scalings {0,0,0,1,0} 
and {0,1,0,0,0}. The same happens in loop-momentum 
space, where the NLO Glauber contribution is proportional 
to (A3 — A2XA5 — A4), thus vanishing in the case considered 
here. So the NLO contribution to F(Q 2 ,m 2 ) is made up en- 
tirely from the LO collinear contributions. The 1 -collinear 
region provides 

_. d/2 e -i*{h+2h+2h) e j^r{X 2 - A 4 )r(i - 2A 4 ) 
~ lK 2(Q 2 ) 2 + 2 ^(m 2 y+ e +^+ 2 hr(i+Xi)r(i+x 2 ) 

r{l + £ + X l +2X 2 )r(-e-2X 2 ) 1 



r(i-A 4 )r(-£-2A 4 ) 1+2A4 



(31) 



in agreement with the momentum-space expansion. The 2- 
collinear contribution is obtained from this by exchanging 
A2 <-> A 4 . Adding the two collinear contributions together 
and performing the limit Ai,A2,A 4 — s- 0, the result (23) is 
reproduced. 

Considering finally the NNLO contributions to 
F(Q 2 ,m 2 ), we expect NNLO Glauber contributions, NLO 
collinear contributions and a LO hard contribution. But here 
the NLO collinear contributions vanish exactly for A3 = A2 
and A5 = A4, both in loop-momentum space and when 
expanding the parametric integrals. So we are left with the 
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NNLO Glauber contributions yielding 



in 



d/2 ie -i^+^2+2X,) r fa _ l) r fa _ l ) 

16(g 2 ) 3 ( OT 2)£+A 1 +2X 2 +2A 4 

r(e + Ai + 2A 2 + 2A 4 )r(l -e-2A 2 -2A4) 

x r(i +Ai)r(i + A 2 )r(i +A 4 )r(-i - e) 

and the LO hard contribution, 

J - e -iK(A 1 +2A 2 +2^t) e i7c(£+A I +2A 2 +2A 4 )/2 



(32) 



/7T 



V2. 



20r(2e 2 ) 3 + £ + ;i i+ 2;i 2+ 2;i 4r(l +Aa)r(l + A 4 ) 

p ^ -l-£-Ai-2A 2 ^ p ^ -l-£-A|-2A 4 

X r(- 1 - 2e - Ai - 2A 2 - 2A 4 ) 

x r ^ -2-£-A)-2A 2 -2A 4 ^j p ^ 3+e+Ai+2A 2 +2A 4 ^ ^ 

both consistent between the expansions of the loop integral 
and of the parametric integrals. Setting all A, = in (32) 
and (33), we recover the results from (24) and (25), respec- 
tively. 

We may also evaluate the parametric integral (29) in 
terms of a onefold Mellin-Barnes representation: 



-lit 



d/2 _ 



-rt(A!+2A 2 +2A 4 ) 



47rr(l+A 1 )r(l+A 2 )r(l+A 4 ) 
2%i J 

r(l + Ai + z )r(3 + e + Ai + 2A 2 + 2A 4 +z) 



r(-l - 2e - Ai - 2A 2 - 2A 4 - z) 



x r(-z)r 



-l-£-Ai-2A 2 -c 



— 1— £— A[— 2A 4 — z 



xr 



2 / -2-£-A|-2A 2 -2A 4 -; 



(34) 



The relevant regions can easily be determined from the 
gamma functions in the last two lines. In particular, the 
squared gamma function indicates that the expansion of the 
parametric integrals I± requires two regions for the Glauber 
contribution, both scaling as (m 2 )~ 2 ~ £ ~ Xl ~ 2Xl ~ 2X4 at lead- 
ing order. When combining F(Q 2 ,m 2 ) = 2{I + +/_), one of 
these gamma functions is cancelled, and we obtain 



F(Q 2 ,m 2 ) = in d/2 



,- e -i7r(A 1 +2A 2 +2A 4 ) e in:(£+A 1 +2A 2 +2A 4 )/2 

20rr(l +Ai)r(l +A 2 )r(l +A 4 ) 



x — f dz(m 2 Y(2Q 2 )- 3 - £ - Xi - n 2- 2X i- z e iKZ / 2 
2ni J 



r(i+A!+z)r 



3+£+A|+2A 2 +2A 4 +^ 
2 



r(-l -2e-Ai -2A 2 -2A 4 -z) 



xr(-z)r 



-l-e-Ai-2A 2 -; 



r 



-l-£-Ai-2A 4 -- 



Xf 



-2-£-A]-2A 2 -2A 4 -z 



(35) 



From this representation the contributions to the asymptotic 
expansion in the limit m 2 /Q 2 — > can be extracted: The hard 



contributions stem from the residues of the poles at z = «, 
the 1- and 2-collinear contributions from z = — 1+2« — e — 
Ai — 2A 2 and z = — 1 + 2n — £ — Ai — 2A 4 , respectively, the 
Glauber contributions from z = — 2 + 2n — e — Ai — 2A 2 — 
2A 4 (with n = 0,1,2,.. .). All LO, NLO and NNLO contri- 
butions reported in (30)-(33) are confirmed by the corre- 
sponding residue contributions from (35). In particular, the 
structure of the poles in (35) clearly shows that the Glauber 
region does not contribute to the NLO result, and that the 
collinear contributions are absent at NNLO, as obtained be- 
fore. 

The results reported in this section show that the use of 
generic propagator powers helps disentangling the Glauber 
and collinear regions from each other, making all regions 
contribute in the same way to the asymptotic expansion in 
loop-momentum space and to the expansion of the paramet- 
ric integrals. Keeping the dependence of the contributions 
and their scalings on the propagator powers also facilitates 
the identification of regions found by asy2 . m in paramet- 
ric space for a subsequent expansion at the level of the loop 
integration. 



6 Summary of asy2 . m 

The updated version of the code, asy2.m, can be down- 
loaded from the known web site [18], where further in- 
stallation instructions are found. The Mathematicacode is 
loaded using <<asy2 .m. 

The main function AlphaRepExpand [] identifies all re- 
gions which contribute to the asymptotic expansion of a 
given loop integral: 

AlphaRepExpand [{kl , k2 , . . . } , 

{(kl + pl)~2 - ml~2, (k2 + p2)~2 - m2~2, 
...}, 

{pl~2 -> Ql, p2~2 -> Q2, pl*p2 -> Q3, ...}, 
-[ml~2 -> x, m2~2 -> x~2, 

Ql -> 1, Q2 -> 3/2, . . .}, 
options] 

The first argument is the list of loop momenta. The sec- 
ond argument lists the denominators of the propagators. The 
third argument contains replacement rules for all kinematic 
invariants. In particular, all external momenta appearing in 
the denominators must be replaced here, otherwise they are 
not identified correctly as vectors. The fourth argument sets 
the scaling of the parameters by replacing all symbols with 
powers of the expansion parameter, labelled by the global 
symbol x, and rational numbers (integers or explicit frac- 
tions of integers). 

The output is a list of regions, specified by the scaling 
(in powers of the expansion parameter) of the Feynman pa- 
rameters jci,JC2, ... corresponding to the propagators in the 
order stated in the second argument. E.g. the output 
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{..., {0, 2, 1, ...}, ...} 

indicates that there is a region {0,2, 1, . . .} specified by the 
scaling x\ ~ x°, X2 ~ x 2 , X3 ~ jt 1 , ... of the Feynman param- 
eters, where x is the small parameter of the problem. 
Possible options of AlphaRepExpand [] are: 

- PreResolve -> True: Try to eliminate cancellations 
between terms in the parametric representation by de- 
composing the integral and performing changes of vari- 
ables. The output contains for each region a list of three 
entries, i.e. a region is e.g. specified by: 

«x[l] -> y[l]/2, x[2] -> y[l]/2 + y[2]}, 
2, {0, 1/2}} 

The first entry in this list is the transformation between 
the original Feynman parameters x[l], x[2], ... and 

the new variables y [1] , y [2] , The second entry is 

the Jacobian of the integral transformation (here "2", i.e. 
dyidj2 = 2ckick2). The last entry specifies the region in 
the usual form with the scalings of the new variables in 
powers of the small parameter, here y\ ~ x°, y2 ~ x 1 / 2 . 
When the option PreResolve is enabled, the code 
warns if it fails to eliminate all possible cancellations 
in the parametric representation. When no warning is is- 
sued, all regions are found. Without this option, how- 
ever, regions occurring at such cancellations will not be 
revealed. 

- GenericPowers -> True: Take into account generic 
(in particular non-integer) powers of the propagators, 
e.g. when these powers are used as analytic regulators. 
Without this option, the preresolution algorithm trig- 
gered by the option PreResolve only finds all regions 
for integrals with propagators raised to positive integer 
powers. 

- Verbose -> True: Print verbose internal information. 

- Scalar -> True: Permit more complex structures of 
the denominators by specifying scalar products of mo- 
menta via the function Scalar [k,p] instead of simple 
products k*p or Scalar2 [k] instead of k~2, e.g.: 

AlphaRepExpand [{k} , 
{Scalar [k, k] - m~2, 

Scalar2[k - q] - m~2}, 
{Scalar2[q] -> qq, m~2 -> qq/4 + y}, 
{qq -> 1, y -> x}, 

PreResolve -> True, Scalar -> True] 

For expanding more general integrals, which need 
not originate from Feynman diagrams, the command 
WilsonExpand [] may be used: 

WilsonExpand[F, U, {xl, x2, ...}, 
{... -> x, ...}, options] 



Traditionally, the first two arguments are the polynomials 
T and U, respectively, from the alpha parametric repre- 
sentation (3) of the Feynman integral or from the Feyn- 
man parametric representation (4), (5). But, more generally, 
WilsonExpand [] reveals regions for integrals over param- 
eters X\,X2, ■ ■ ■, integrated from to °° each, where all non- 
trivial polynomials of the parameters x, occurring in the inte- 
grand are specified either in the first or the second argument. 
The third argument of WilsonExpand [] is the list of inte- 
gration parameters. The fourth argument specifies the scal- 
ing of all quantities with the expansion parameter x, in the 
same way as in the fourth argument of AlphaRepExpand [] . 

The output is a list of regions, specified by the scal- 
ing of the integration parameters in powers of the ex- 
pansion parameter, exactly as described for the output of 
AlphaRepExpand [] (without the PreResolve option). 

Possible options of WilsonExpand [] are: 

- Delta -> True: Under the integral, the sum over an 
arbitrary non-empty subset of the integration parameters 
is restricted to 1 via a delta function. The specific choice 
of the sum in this delta function must be irrelevant for the 
integral, which is the case for the generalized Feynman 
parametric representation (4). 

- Verbose -> True: Print verbose internal information. 

More generally, the option Delta of WilsonExpand [] 
works correctly for all integrals of the form 

J o ■■■J o dxi---dx N S (^Y/^iXi-lj /(xi,...,xn) , (36) 

where the linear combination in the argument of the delta 
function has no negative and at least one positive coefficient 
(a, > 0V7 and 3 a, > 0), and where the function / scales ho- 
mogeneously with the parameters Xi as 

f{Xx\ , Xxn) = X~ N f{x\ ,xn) VA > , (37) 

the degree of homogeneity being equal to minus the num- 
ber of integration parameters. It can be shown that for such 
integrals (36) the specific choice of the coefficients a, is ir- 
relevant. 5 The integrand (5) of the Feynman parametric rep- 
resentation (4) and the integrands of all parametric repre- 
sentations of Feynman integrals used in this paper (without 
the additional analytic regularization factors jc?') fulfill the 
homogeneity condition (37), so they do not depend on the 
specific choice for the arguments of their delta functions. 

To see this, multiply the integrand of (36) by 1 in the form 
Xj JjJ° dt e~ ,x i , where Xj is any of the integration parameters. Then, 
inside the f-integration, transform the integration variables as 
Xi — > Xj/t, i = l,...,iV, and use the homogeneity relation (37) 
with X = 1/f. Finally evaluate the f-integration first, yielding 
Jq dt S 1 a i x i ~t) = 1) independent of the coefficients as long 
as the lineal' combination is positive. The integral (36) is given by 
/o°° ' ' ' So d*i ' ' ' f{x\ , ■ ■ ■ , x N ) Xj e~ x i . 
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7 Conclusion 

We have presented an algorithm for identifying all regions 
which are relevant for the asymptotic expansion of a given 
loop integral at the level of its parametric representation. In 
contrast to the previous version asy . m of the code, also po- 
tential regions and Glauber regions are found now. The nec- 
essary decompositions and variable transformations of the 
integral are automated by the updated Mathematica code 
asy2.m [18]. When the command AlphaRepExpand [] is 
used with its option PreResolve enabled, we are convinced 
that it either reveals all relevant regions or issues a warn- 
ing. In particular, regions corresponding to cancellations be- 
tween large positive and negative terms in the parametric 
representation of the loop integrals (such as potential and 
Glauber regions) will now be found. 

Let us emphasize that to prove expansion by regions at 
least for some specific limit typical of Minkowski space is a 
natural mathematical problem. Perhaps, this problem is not 
specifically related to Feynman integrals. Let us present an 
example of a one-dimensional parametric integral, without 
any relevance to Feynman integrals, and show that expan- 
sion by regions works successfully. To do this, we will use 
asy2 .m. 

Let us consider the integral 



F (t)= f (t + u + u 2 ) l du 
Jo 



(38) 



with A a complex parameter, in the limit t — > 0. We assume 
that A is in the domain ReA < — 1/2 in order to have an ab- 
solute convergence of the integral which then can be contin- 
ued analytically to the whole complex plane as an analytic 
function of A . Running 

WilsonExpand[t + u + u~2, 1, {u}, {t -> x}] 

we obtain the two regions {{1}, {0}}. The leading-order 
terms from each region can be evaluated analytically in 
terms of gamma functions at general A, with the results 

f A +'r(-A- 



r(-A) 

and 

r(-2A-i)r(A + i) 
tFa! 

They can be checked easily by deriving the onefold Mellin 
Barnes representation 

1 1 

F(t) 



(39) 



(40) 



27nT(-A) 

x /r(-z)r(A-z+l)r(-2A + 2z-l)/*dz (41) 



and evaluating the first terms of the asymptotic expansion 
in the limit t — > by shifting the contour to the right and 
taking residues at the poles of the two gamma functions in 
the integrand. 
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